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Abstract 

We present the explicit expressions for the resummation of large-XT Sudakov effects in the 
transverse-energy distribution of prompt-photons produced in hadronic collisions, to next- 
to-leading logarithmic (NLL) accuracy. Fragmentation processes do not contribute to the 
Sudakov resummation at NLL level. In Mellin space, the resummed radiative factor factorizes 
in the product of independent radiators for the two initial and the one final coloured partons 
appearing in the Born process, times a simple factor describing the soft-gluon interferences 
between initial and final states. The formulae are given in terms of Mellin moments, and can 
be used for phenomenological applications using standard techniques for the inverse-Mellin 
transforms. The calculations presented in this work, when added to the existing works on 
DY and DIS production, complete the theoretical ground-work necessary to carry out global 
fits of parton densities with a uniform NLL accuracy in the large- x region. 
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1 Introduction 



Prompt-photon production plays a very important role in our understanding of the physics of 
hadron collisions. At the leading order (LO) in QCD perturbation theory, prompt photons are 
produced via light-quark annihilation, with emission of a hard gluon recoiling against the photon, 
or via quark-gluon Compton scattering, with the emission of a quark. When the Compton process 
dominates the cross section, tests or even measurements of the gluon density inside the proton 
can be performed. This is the case of photon production at small xt = 2E T /\fS (xt & 0.1) in 
pp collisions, and it is true for all values of xt that are accessible in fixed-target proton-nucleon 
collisions, due to smaller content of antiquarks relative to gluons in the nucleon sea. In particular, 
prompt-photon production at large Xt can therefore be used to constrain or measure the gluon 
density at large x. This region of the gluon density is of great importance for the study of 
high-transverse-momentum phenomena at hadron colliders, and it is not accessible using only 
Deep-Inelastic-Scattering (DIS) data. 

The cross section for inclusive photon production has been computed at the next-to-leading 
order (NLO) in perturbation theory M— [EJ. The NLO computation of isolated photon produc- 
tion 0-0, which is the relevant quantity for the measurements carried out in high-energy pp 
collisions, is also available in the limit of small size of the isolation cone 0-0. These calculations 
include all light-parton fragmentation processes up to NLO |10|| . 

These theoretical results have been used |]11||-|]14|| to probe the overall consistency of the 



prompt-photon production data from both fixed-target |15|-|16| and collider experiments [JT7 
P 



The interpretation of the data has not provided so far a fully satisfactory picture. The study by 
Huston et al. fT3| exposed a tendency of the xt distributions to be steeper than theory, regardless 
of the value of \/S. This result could not be accomodated by a simple modification of the parton 
densities, since different experiments probe different values of xt- These authors therefore pro- 
posed that additional mechanisms should be introduced to explain the pattern of the data. In the 
fixed-target regime, such a mechanism would be provided by the presence of a non-perturbative 
kx kick, which would smear the Et spectra. This phenomenon was also advocated to help explain 
the spectra of fixed-target heavy-quark production [pOR . In the high-energy regime, probed by the 
Tevatron experiments, the Et smearing necessary to reconcile theory and data could be provided 
by the inclusion of multi-gluon emission effects from the evolution of the initial state, as advocated 
by Baer and Reno |pl|| . 



The analysis by Vogelsang and Vogt [13] indicated that allowing for different choices of fac- 
torization and renormalization scales, a satisfactory fit to the data could be accomodated by 
modifying the gluon density within the range allowed by the DIS data available in 1995. This 



interpretation is apparently not viable anymore p2| , because of the most recent constraints on 
the gluon density extracted at small x from the HERA data, and, in particular, because of the 
latest fixed-target prompt-photon data from E706 



The comparison of the E706 data with NLO QCD, carried out in |23|, seems to confirm the 



need for an intrinsic- kx smearing corresponding to (kx) ~ 1 GeV. These conclusions are shared in 
a recent global fit of the parton densities performed by the MRST group |24[] . In this same study 
(see also p5]), however, a strong dependence of (k T ) on v^S is claimed to be necessary to properly 
describe the lower-energy data published by WA70. 

In conclusion, the comparison of the large-x^ fixed-target data with NLO QCD still presents 
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some puzzling features, which will need to be properly clarified before use of these data can be 
made to place robust constraints on the large- x gluon distribution inside the proton. This is 
unfortunate, since these data provide today the only independent probe on high- a; gluons. Their 
accurate interpretation is therefore a fundamental ingredient for an accurate prediction of the 
production rate of high-i?T jets at the Tevatron, a measurement which has challenged perturbative 
QCD in the recent past f26fl . 

To improve the reliability of the perturbative predictions for the production of prompt photons 
at large xt, and detect the presence of potentially large corrections beyond NLO that could 
change the interpretation of the current data, in the present work we consider an extension of the 
NLO formalism that includes large logarithmically-enhanced effects as the production threshold 
is approached. This is the kinematical region of interest for the fixed-target data. As the xt of 
the photon is increased, the parton luminosity becomes steeper, being driven down by the strong 
suppression of the gluon density at large x. We thus enter a regime of inhibited radiation: further 
radiation of soft gluons is strongly suppressed, and logarithmically-enhanced effects (Sudakov 
effects) arise at any order in the perturbative expansion. These logarithms spoil the reliability of 
the fixed-order expansion in the strong coupling a s and, hence, their summation to all pertubative 
orders is necessary. For simplicity, this regime can be described in terms of the distance from 
the kinematic threshold, which is reached when xt ~ 1. In this limit, the coefficients of the 
perturbative series for the cross section are enhanced by powers of m(l — xt) that have to be 
resummed at all orders. This simplified description applies to the case of hypothetical structure 
functions that are not strongly suppressed at high x. It is, however, an appropriate framework for 
the classification of the perturbative corrections we are interested in. 

In this work we will present all the formalism that is needed to compute the resummed cross 
section for direct photons, integrated over the photon rapidity and at fixed transverse energy. In 
particular, we give explicit resummation formulae that are valid up to next-to-leading logarithmic 
(NLL) accuracy. No phenomenological applications will be discussed here, but they will be ex- 
plored in a forthcoming work. Furthermore, the formulae for the resummed correction factors will 



be presented and illustrated, but not derived here. The general formalism [2/] used to obtain the 



resummation factors, which has already been used for the NLL resummation of the heavy-quark 



total production cross-section in Ref . |5q| , will be presented in a forthcoming publication 

The rest of this work is organised as follows. The general theoretical framework is discussed in 
Sect. 2. In Sect. 3 we fix our notation and present the formulae for the Born cross section, together 
with their Mellin transforms. Soft-gluon resummation at large xt is considered in Sect. 4. The 
NLL resummation factors are presented in Sects. 4.1, 4.2. In Sect. 4.3, the fixed-order expansion of 
the resummed formulae is compared with the NLO results of Refs. |L], [|. This comparison is also 
exploited to fix certain constant factors in the resummed formulae. In Sect. 5 we discuss similarities 
and differences between the resummed factors for the prompt-photon cross section and those for 
other hard-scattering processes, and we prove the consistency of the results obtained in the case 
of prompt photons with the coherence properties of large-angle soft-gluon emission. Section 6 
contains our conclusions. 

More technical details are left to the Appendices. In Appendix A, we give the NLL formulae 
for the radiative factors in the Mellin transform representation. Previous experience in the case of 
heavy-flavour production has shown that this is what is needed to perform a reliable phenomeno- 
logical analysis f28f . In Appendix B, the threshold limit of the partonic cross sections is discussed. 



In particular, a prediction for the logarithmic terms at the next-to-next-to-leading order (NNLO) 
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is given. Finally, in Appendix C the simpler case of photoproduction of direct photons is discussed. 

While completing this paper, a study of the NLL resummation for single-inclusive distribu- 
tions, covering the case of prompt-photon production, has been released by Laenen, Oderda and 
Sterman [jSGfl . 



2 General framework 

The presence of large logarithmically-enhanced contributions is a common feature in the study of 
the production cross sections of systems of high mass or high transverse energy near threshold. 
In this kinematic regime, known as the Sudakov regime, only additional soft gluons can be pro- 
duced. The radiative tail of the real emission is thus strongly suppressed and cannot balance the 
virtual corrections. The imperfect compensation between real and virtual terms leads to the large 
logarithmic contributions. 

General techniques for resumming soft-gluon corrections to hadropro duct ion processes have 
been developed over the past several years, starting from the case of Drell-Yan (DY) pair pro- 
duction |3l], |32[ . The resummation program of the soft-gluon contributions is best carried out in 
the Mellin-transform space, or iV-space, where N denotes the parameter that is conjugate to the 
kinematic variable that measures the distance from threshold. In iV-moment space the threshold- 
production region corresponds to the limit iV — > oo and the typical structure of the logarithmic 
contributions is as follows 



a 



(0) 
N 



2)i 



n=l m=l ) 

onic cross-section at LO. 
terms in the curly bracket of Eq. (ffl) can be explicitly summed and organized in a radiative factor 



where a$ is the corresponding partonic cross-section at LO. In the DY process the logarithmic 



A dy,n that has an exponential form |3T|, |32|, |33| : 

{oo n+1 ^ 

J2<J2 G nin^ m N\ (2) 
n=l m=l ) 

= exp{lnNg i i Y {a s \nN)+g§ Y (a s lnN)+a s g^ Y {a s lnN)+...} . (3) 

Note that the exponentiation in Eq. @ is not trivial. The sum over m in Eq. (P extends up to 
m = 2n while in Eq. (0) the maximum value for m is smaller, m < n + 1. In particular, this 
means that all the double logarithmic (DL) terms OL n s c n ^ n ln 2n N in eq. ([!]) are taken into account 
by simply exponentiating the lowest-order contribution ol s C\^ In 2 N. Then, the exponentiation in 
Eq. (0) allows one to define the improved perturbative expansion in Eq. ([3]) . The function In N g^ Y 
resums all the leading logarithmic (LL) contributions a" ln n+1 N, g^ Y contains the next-to-leading 
logarithmic (NLL) terms a"ln n N, Oi s g DY contains the next-to-next-to-leadingloga.iith.mic (NNLL) 
terms ln n N, and so forth. Once the functions g^ Y have been computed, one has a systematic 
perturbative treatment of the region of N in which a s lnA^ 1, which is much larger than the 
domain a s In 2 iV 1 in which the fixed-order calculation in a s is reliable. 

The QCD exponentiation formula for the DY process formally resembles analogous results for 
QED. This is because the underlying hard-scattering subprocess involves only two QCD partons, 
i.e. the annihilating qq pair, and, hence, both its kinematics and its colour structure are simple. 
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In the case of prompt-photon production, instead, all the LO hard-scattering subprocesses 



q + q -> g + 7 , q + g q + i , q + g ^ q + i , (4) 

involve three coloured partons and, then, Sudakov resummation is by far less trivial. 

A key ingredient for the exponentiation in the DY process is the factorization of the corre- 
sponding multigluon matrix elements in the soft limit. Since the colour structure of a two-parton 
hard-scattering is trivial^], primary soft radiation from the two hard partons factorizes as in QED. 
Then the subsequent parton radiation can be factorized in non-interfering angular-ordered cas- 
cades because of the coherence properties |35[ of QCD emission. 



In the case of prompt-photon production, and, in general, in scattering processes produced by 
hard interactions of more than two QCD partons, the colour and momentum flows in the partonic 
subprocess are more involved. In particular, the interplay between colour exchange in the hard 
scattering and colour transitions induced by parton radiation spoils QED-like factorization of 
soft-gluon emission. Therefore, both colour correlations and soft-gluon interferences have to be 
properly taken into account. It turns out that, in general, the threshold logarithmic corrections 



cannot be resummed in a single exponential factor [ : one has to deal with exponential matrices 
that couple the various colour channels of the hard-scattering subprocess. 

However, the three-parton subprocesses in Eq. (|J) are a special case among the multiparton 
configurations. There is only one colour-singlet stateJJ] that can be constructed by combining qqg 
and then, because of colour conservation, soft-gluon radiation cannot induce colour transitions in 
the hard-scattering subprocess. Owing to the absence of colour correlations, we conclude that the 
logarithmically-enhanced threshold corrections in prompt-photon hadroproduction are embodied 
by three radiative factors (one factor for each of the LO partonic channels in Eq. (^)) 

A*T ff > s ), AST^a.), A?T 9 "> S ), (5) 

that, after all-order resummation, have an exponential form analogous to the DY radiative factor 
in Eq. (||). Nonetheless, the similarity with the DY process regards only the colour structure. The 
hard-scattering kinematics is different in prompt-photon production and the factors in Eq. (|5|) 
still contain soft-gluon interference effects that are non-trivial. The pattern of these soft-gluon 
interferences is typical of multiparton hard-scatterings. 

General theoretical methods to perform Sudakov resummation in processes initiated by hard 
scattering of more than two QCD partons have recently been developed by two groups. The 
KOS formalism ||36||-||38|| uses the Wilson line approach to treat colour correlations and soft-gluon 
interferences. It has been explicitly applied to the calculation with NLL accuracy of the the 



invariant-mass distributions of heavy-quark pairs and dijets (see also ref. |[30|| ). The more recent 



BCMN formalism |28| is based on generalized soft-gluon factorization and has been used for the 
NLL calculation of the total cross section for heavy-quark hadroproduction. The consistency of 
the NLL results for the total cross section with those for the invariant mass distribution |36] 



of heavy-quark pairs shows that, although different, the two formalisms are equivalent to a large 
extent. 



^This is the reason why similar exponentiation formulae apply to many other <wo-je<-dominated processes |34| . 

"in other words, the qqg colour-amplitude -M^g (ce, a and a are the colour indices of the quark, antiquark and 
gluon, respectively) is necessarily proportional to the matrix of the fundamental representation of the gauge 
group .S7 •: A.. :. 
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In the rest of this paper, we first introduce our notation and then we present the resummed 
expressions of the prompt-photon radiative factors (|5]) to NLL accuracy. The results include the 
complete soft-gluon interferences to this accuracy, as evaluated by using the BCMN formalism. 



Details of our general formalism will be presented elsewhere [29 



3 Notation and fixed-order calculations 

We consider the inclusive production of a single prompt photon in hadron collisions: 

H X {P X ) + H 2 (P 2 ) - 7 (p) + X . (6) 

The colliding hadrons Hi and H 2 respectively carry momenta P" and P 2 . In their centre-of-mass 
frame, using massless kinematics, they have the following light-cone coordinates 

i? = yf (1,0,0), ^ = ^1(0,0,1) , (7) 

where S = {P\ + P 2 ) 2 is the centre-of-mass energy squared. The photon momentum p is thus 
parametrized as 

where Et and r\ are the transverse energy and the pseudorapidity, respectively. We also introduce 
the customary scaling variable xt (0 < xt < 1): 

We are interested in the prompt-photon production cross section integrated over rj at fixed Et- 
According to perturbative QCD, the cross section is given by the following factorization formula 



<icr 7 (xT, E T ) 

dE T rir „.,, 



1 f 1 f 1 

—3-^/ dx x fa/H^Xi, (1%) / dx 2 f b/H2 (x 2 ,H 2 F ) 

[ dxldix ^== ] dr ab ^(x, a s (fi 2 );ET,fi 2 ,fi 2 F ,fi 2 f ) (10) 

JO [ \ y/XlX 2 J J 

+ V/ dzZ 2 d c /Jz,lA)8[x ^T__ \ J^^^^).^ 2^2_ 2j 1 _ 

c J J V Z^/XiX 2 ) 3 J 

where a,b,c denotes the parton indices (a = q,q,g), and f a /Hi( x i, a40 an d fb/H 2 ( x ij are the 
parton densities of the colliding hadrons evaluated at the factorization scale /ip. The first and 
the second term in the curly bracket on the right-hand side of Eq. ( |l0l) respectively represent 
the direct and the fragmentation component of the cross section. The fragmentation component 
involves the parton fragmentation function d c / 7 (z, fi 2 ) of the observed photon at the factorization 
scale fif, which, in general, differs from the scale \xf of the parton densities. 
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The rescaled^ partonic cross sections a ab ^ 7 and a a b->c in Eq. ( PH) are computable in QCD 
perturbation theory as power series expansions in the running coupling a s (fi 2 ), \x being the renor- 
malization scale in the MS renormalization scheme: 



a ab ^ c (x, a s (fi 2 ); E T , [i 2 , /i F , f/ f ) = a 2 (fi 



°abLd-y( X ) + a W) °abLy& e t , /i 2 , /4, n)) 

n=l 

(ID 

oo 

2\. iti2 ,,2 ,,2 ,.2\ „2/..2\ 



(12) 



°abLc(x) + «r(/^ 2 ) VabLcfa E T, /A VF, Vf) 
n=l 

Note that the ratio between the direct and the fragmentation terms in Eqs. (|TT]) and (|l|) is of the 
order of a/a s , where a is the fine structure constant. This ratio is compensated by the photon- 
fragmentation function eL/ 7 , which (at least formally) is of the order of a/a s , so that direct and 
fragmentation components equally contribute to Eq. (]10D . 

Throughout the paper we always use parton densities and parton fragmentation functions as 
defined in the MS factorization scheme. In general, we consider different values for the renormal- 
ization and factorization scales fi,fip,fif, although we always assume that all of them are of the 
order of the photon transverse energy Et- 

The LO terms cr^ d7 in Eq. ([□]) are due to the tree-level parton scatterings 

a + b-> + 7 , (13) 

where the flavour indices a, b, d are those explicitly denoted in the subprocesses of Eq. (|J). Using 
our normalization, the two independent (non-vanishing) partonic cross sections for the direct 
component are: 



1 X 



2 



*£L» = *2U(*) = *<W C 7fr=p [ l + tJ ' (15) 

where e q is the quark electric charge. Note that, having integrated over the photon pseudorapidity, 
the expressions (|14|, 15) are even functions of the photon transverse energy Ex, i.e. they depend 
on x 2 rather than on x. The NLO terms m Eq. (pT|) were first computed in Ref. flj. 

The partonic contributions cfab-+c to the fragmentation component of the cross section are 
exactly equal to those of the single-hadron inclusive distribution. Their explicit calculation up to 
NLO was performed in Ref. [|TIJ . 

We are mainly interested in the behaviour of QCD corrections near the partonic-threshold 
region x — ► 1, i.e. when the transverse energy E? of the photon approaches the partonic centre- 
of-mass energy ^/xix^S. In this region, the LO cross sections ([14], Il5|) behave as 



*2U(*) - 7fTT2 ■ (is) 

This integrable singularity is a typical phase-space effect. At higher perturbative orders, the 
singularity in Eq. ([16]) is enhanced by double-logarithmic corrections due to soft-gluon radiation 
and the cross section contributions in Eqs. ( |TT] , |T^ ) behave as 

aW(x)~a®{x) ln 2n (l-x) . (17) 



*Thesc functions are related to the partonic differential cross sections by cr a b^i — E\da a b^i I AEt {% = 7, c) 
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Resummation of these soft-gluon effects to all orders in perturbation theory can be important to 
improve the reliability of the QCD predictions. 



3.1 N- moment space 



The resummation program of soft-gluon contributions has to be carried out [|31], |32| in the Mellin- 
transform space, or iV-space. Working in iV-space, one can disentangle the soft-gluon effects in the 
parton densities from those in the partonic cross section and one can straightforwardly implement 
and factorize the kinematic constraints of energy and longitudinal-momentum conservation. 

It is convenient to consider the Mellin transform <j 1 ^{Et) of the dimensionless hadronic distri- 
bution E^dcT 1 {xTi Et)/cIEt- The iV-moments with respect to x\ and at fixed Et are thus defined 
as follows: 

a.MEr) = T dxl (4)^ j ^M^r) (lg) 
Jo clEt 

In iV-moment space, Eq. (|T0|) takes a simple factorized form 



^,n{Et) — ^2 fa/H u N+l(^F) fb/H 2 ,N+l{l^F) 
a.b 

• {<5-a6^ 7 , Af(«s(/i 2 );-E|,/i 2 ,/Lt|,/i/) (19) 
+ °ab->c, N(a S (l^ 2 )] E^, jj 2 , jJ F , jjj) 4/ 7 ,2AM-3(/4) j > 

where we have introduced the customary iV-moments f a /H,N and d a /7,jv of the parton densities 
and parton fragmentation functions: 

fa/H,N{^ 2 ) = f dx X N ~ X f a/H {x,II 2 ) , (20) 
•/ 

da/ 7 ,N(^ 2 ) = dz z N ~~ 1 d a/j (z, yU 2 ) . (21) 



Note that the iV-moments of the partonic cross sections in Eq. (19) are again defined with 



respect to x\: 



a a6 ^ 7jAr (a s (/i 2 );£|,/i 2 ,/4>/4) = ^ dx 2 (x 2 ) N 1 a a6 ^ 7 (x,a s (/i 2 );£|,/i 2 ,/4,/iJ) . (22) 

In particular, the iV-moments of the LO contributions in Eqs. (|14| , |15| ) are given by the following 
explicit expressions: 

*v> -rr* CF r(i/2)r(iv + i) 

a aq-+q^ N = T B — (2 + N) , (23) 



*(o) _ *<o) _ „ a 2 J_ r(l/2) T(N + 1) 

8N C T(iV + 5/2) 



*2Ur. n = frgL*, n = - e 2 ^ ^ i^Z ^ (7 + 5iV) . (24) 



Note also the pattern of moment indices in the various factors of Eq. (IS), i.e. f a /H,N+\ f° r the 
parton densities and <i c / 7i 27v+3 for the parton fragmentation functions. This non-trivial pattern 
follows from the conservation of the longitudinal and transverse momenta. 
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The threshold region iy — > 1 corresponds to the limit iV — > oo in iV-moment space. In this 
limit, the soft-gluon corrections (\Tf\) to the higher-order contributions of the partonic cross sections 
become 

In 2 " iV . (25) 

The resummation of the soft-gluon logarithmic corrections to all orders in perturbation theory is 
considered in the following Section. 



4 Soft-gluon resummation at high Et 

4.1 Resummed cross section to NLL accuracy 

In the threshold or large- N limit, the various partonic channels contribute in different ways to the 
prompt-photon cross section <t 7 ^ n (E t ) in Eq. (|19|). 

Firstly, we can compare the direct and fragmentation contributions to Eq. ([19]). The partonic 
cross sections 6" a ; ) _ +7i jy and a a b^c, n have the same large- N behaviour, but, owing to the hard 
(although collinear) emission always involved in any splitting process c — > 7 + X, the photon- 
fragmentation function d c / 7j 7v is of the order of 1/N. Therefore, in Eq. flT9|) the fragmentation 
component is formally suppressed by a factor of 1/N with respect to the direct component and in 
our resummed calculation we can neglect the fragmentation contributions. 

Then, we can discuss the differences in the large-iV behaviour of the partonic cross sec- 
tions <r a f,^ 7i Ar(a s ) for the direct processes. The cross sections for the partonic channels ah = 
qq', qq', qq, qq', qq, qq' (q and q' denote quarks of different flavours) vanish at LO and are hence 
suppressed by a factor of a s with respect to (T g ^ 7j 7v(a s ), o qg ^^ Ar(a s ), <3"g g ^ 7 , 7v(a s ). Moreover, in 
the large- N limit this relative suppression is furtherly enhanced by a factor of 0(1/N) because 
the photon has to be accompanied by (at least) two final-state fermions that are not produced by 
the decay of an off-shell gluon. Therefore, we make no attempt to resum soft-gluon corrections 
to these partonic channels. The partonic cross section a gg ^ 1 ^{a s ) has a different large-iV be- 
haviour. It begins to contribute at NLO via the partonic process g + g^'y + q + q, which again 
leads to a suppression effect of 0{1/N) with respect to the LO subprocesses. However, owing to 
the photon-gluon coupling through a fermion box, the partonic subprocess g + g^j + g is also 
permitted. This subprocess is logarithmically-enhanced by multiple soft-gluon radiation in the 
final state, but it starts to contribute only at NNLO in perturbation theory. It follows that the 
partonic channel ah = gg is suppressed by a factor of a 2 with respect to the LO partonic channels 
ah = qq, qg, qg and it enters the resummed cross section only at NNLL accuracy. 

In conclusion, since we are interested in explicitly perform soft-gluon resummation up to NLL 
order, we can limit ourselves to considering the partonic cross sections o" g g^ 7 , o qg ^, tiqg^-y 

As discussed in Sect. |2], the soft-gluon corrections to the partonic channels ah = qq, qg, qg are 
not affected by colour correlations. Thus, in the resummed expressions <3"^!!^ 7 N for the partonic 
cross sections, the logarithmically-enhanced threshold contributions can be factorized with respect 



to the corresponding LO cross sections cr^_^ dl N in Eqs. (|3|, ^4|). The all-order resummation 



formulae are 

£9?%, N( a s(v 2 ); El, fi 2 , /i 2 f ) = a a s (fi 2 ) o-JJL^, 7v C q q^(a s (fi 2 ), Q 2 / n 2 ; Q 2 //4) 

■ A^/> s (/i 2 ), QV/^QV/i 2 ,) , (26) 
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a^ g e ^ N {a s {fi 2 )]E^,fi 2 ,fi 2 F ,fi 2 f ) = aa s {fi 2 ) ffJjL^, n C qg ^{a s {fi 2 ),Q 2 / fi 2 ;Q 2 / fi 2 F ) 

■ a^>.(/AW;Q7/4) , (27) 

y N {a s {fi 2 )]E T ,fi 2 ,fi 2 F ,fi 2 f ) = a ( g g e ^ N (a s {fi 2 )]E T ,fi 2 ,fi 2 F ,fi 2 f ) , (28) 



a 



where 



Q 2 = 2E T . (29) 



Note that the right-hand side of Eqs. (|26| , p7|) does not depend on the factorization scale /i/ 
of the photon fragmentation functions. Thus, the resummed partonic cross sections <3"^!!^ 7 N turn 
out to be independent of ///. This is in agreement with the subdominance of the fragmentation 
contributions near threshold, as discussed above. 

The functions C a fe_» 7 (a s ) in Eqs. ( |26| , |2"T| ) do not depend on N and, thus, contain all the 
contributions that are constant in the large- N limit. These functions are computable as power 
series expansions in a s 



C^ 7 (« s (/i 2 ),Q 2 //^ 2 ;Q 2 /^) = i + E P^ 1 C%l ( (Q 2 /^;Q 2 /ti) . (30) 



+ 0O 

E 

n=l 

The physical origin and the structure of the constant factors C a b^{a s ) is discussed in Sect. |Q 




The In iV-dependence of the resummed cross sections is entirely embodied by the radiative 
factors A^^ dl on the right-hand side of Eqs. ( p6| , p7|) . Note, the mismatch between the moment 



index of the radiative factor and that of 0"ibLd 7 , n'- the former depends on JV + 1, like the parton 
densities in Eq. (|19|). The explicit expressions of the radiative factors are given in the following 
subsection. 



4.2 The radiative factors 

The soft-gluon factors A^^ 7 depend on the flavour of the QCD partons a, b, d involved in the LO 
hard-scattering subprocess a + b — > d + 7. According to the discussion of Sect. ^|, the resummed 
expressions for A a ^ dl have an exponential form. To explain the exponentiation structure and to 
facilitate the comparison with other hadropro duct ion processes, we use a notation similar to that 
in Ref. p2[ and we write the prompt-photon radiative factors as follows 



A a ^(a s ^ 2 ),Q 2 /^Q 2 /v 2 F ) = A a N (a s (fi 2 ), Q 2 /^ 2 ] Q 2 /^ 2 F ) A b N (^ 2 ),Q 2 /^ 2 ;Q 2 /^ 



■ 4(a s (/i 2 ),Q 2 /^)A!; nt) " 7 K(/i 2 ),Q 2 /^) • (31) 

The resummed formulae to NLL accuracy for the various contributions on the right-hand side of 
this equation are presented below. 

Each term A^(a s (/x 2 ), Q 2 /fi 2 ; Q 2 /fi F ) depends on the flavour a of a single parton, on the 
factorization scheme of the parton density f a /H,N(^ 2 F) an d on the factorization scale //p. In the 
MS scheme, we have 



A^(a s ( / i 2 ),g 2 / / i 2 ;g 2 //i|)=exp / dz — / 4r A a (a s (q 2 )) + 0(a s (a s \nN) k )\ 

[Jo 1 - z J/4 q 2 J 

(32) 
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where A a (a s ) are perturbative functions 

A a K) = ^w + (-) 2 4 2) + o( 

7T \ 7T / 

The lower-order terms and are 

^i 1 ^ = 5 = — C a K , 



a' 



(33) 



(34) 



where C a = Cf if a = q, q and C a = Ca if a = g, while the coefficient K is the same both for 
quarks [|3£| and for gluons j|0[ and it is given byQ 



A I 18 6 / 9 ' 



(35) 



The term jfj(a s (n 2 ),Q 2 / fi 2 ) depends on the parton flavour d and is independent both of the 
factorization scale and of the factorization scheme: 



4(a s (//),Q7^ 



Z N-1 _ 

exp < / dz 

Jo 1 — z 



+ 0(a s (a s \nN) k ) 



1 



>2 2 ^A d (a s (q 2 )) + -B d (a s ((l-z)Q 2 )) 
(i-z) 2 q 2 q 2 2 



(36) 



The functions A d (a s ) are given in Eq. (^) and the functions Bd(a s ) have analogous perturbative 
expansions: 

(37) 



B d {a s ) = + 0(a s 2 ) 



with 3S, 40 



R (i) _ 3 r 
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B d % = --{llC A -AT R N f) 



1 



(38) 



Likewise Jjy, the remaining contribution A^ r nt ' ) in Eq. fl3~l|) is independent of the factorization 
scale and scheme. Nonetheless, it depends on the flavours of all the QCD partons entering the 
LO scattering subprocess: 



a (int) ab— ><i7 
N 



( -i ,JV-1 _ i 

(a s (fi 2 ),Q 2 /fi 2 ) = exp lj o dz — A^, 7 (a s ((l - ^) 2 Q 2 )) + (9(a s (« s In iV) fc ) 



The function .D^^^^s) has the following perturbative expansion 



7T 



with 



D 



(i) 



(C a + C 6 -C d ) In 2 



(39) 
(40) 
(41) 



The factorized structure in Eq. (|31|) has a direct physical interpretation. The factors 
and A^r take into account soft-gluon radiation emitted collinearly to the initial-state partons. 
Consistently, these are the sole factors that depend on the factorization scale fip of the parton 

ttln SU(N C ) QCD, the colour factors are C F = (N^ - 1)/(2N C ),C A = N c and T R = 1/2. 
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densities of the colliding hadrons. The factor Jfj is due to collinear (either soft or hard) radiation 
in the final-state jet that is produced by the fragmentation of the parton d recoiling against the 
triggered photon. The factor A^ nt ' ) contains the contribution of soft-gluon emission at large angle 
with respect to the direction of the hard partons entering the LO scattering subprocess. This 
factor thus embodies the soft-gluon interference effects anticipated in Sect. £| 

According to this interpretation, the perturbative functions in Eqs. ([33], |37|, [40]) measure the 



intensity of the coupling of i) soft- collinear gfuons (function A a (a s )), ii) hard- collinear partons 
(function B a (a s )) and Hi) large-angle soft gluons (function Z) a 6^d 7 (a; s )). Note that, due to their 
collinear nature, the functions A a (a s ) and B a (a s ) depend on the colour and flavour of the sole 
parton a. On the contrary, ^^(as) depends on the colour charges of all the QCD partons. 

The physical origin of the several contributions on the right-hand side of Eq. fl3l]) is furtherly 
discussed in Sect. where we compare the promp-photon radiative factors with the analogous 
resummed factors that control the threshold behaviour of other hadroproduction processes. In 
the rest of this section we limit ourselves to comment on few additional features of the resummed 
contributions to the prompt-photon cross section. 

The various factors in Eq. (|3~TD contribute to the resummed prompt-photon cross section at 
different level of logarithmic accuracy. If we simply consider the double-logarithmic (DL) approx- 
imation, which consists in resumming only the terms a™ ln 2n N, we can neglect the interference 
factor and the B(a s ) function in Eq. (|36|) and we can expand the exponent in A% and Jfj 

to its first order in a 



3 ■ 



A%(a s ( f i 2 ),Q 2 / f i 2 ;Q 2 /n 2 F ) ~ exp {+ 2C a ^- In 2 iv} , (42) 
J d N (a s ^ 2 ),Q 2 /fi 2 ) ~ expj-0^1n 2 iv} . (43) 



The complete set of LL terms is obtained by neglecting the functions B(a s ), D(a s ) in Eqs. fl36[ , 
[39^), by truncating A a (a s ) to their first order and using the LO running of the coupling ct s (q 2 ). At 
the NLL order, also the contribution of the coefficients A^f 1 , B$ and D^_^ d has to be included. 

Note that different scales, e.g. q 2 , (1 — z) 2 Q 2 , (1 — z)Q 2 , appear on the right-hand sides of 
Eqs. ([321 , |36| , |39|) . In particular, the scales in the g 2 -integration limits of Eq. (|32"D are different 
from those of Eq. (^), and the B function in Eq. (|36| ) depends on a s ((l — z)Q 2 ) while the D 
function in the interference contribution fl3~9|) depends on a s ((l — z) 2 Q 2 ). These scales follows 
from the hard-scattering kinematics, which affect in a different way initial- or final-state emission 
and collinear or soft radiation. 

Note, also, that the renormalization scale fi does not explicitly enter the right-hand side of 



Eqs. (|32| , |36| , p9|). This is because the radiative factors are renormalization-group-invariant quan- 
tities when evaluated to all order in perturbation theory. Only when the all-order expressions 
are truncated to a certain degree of logarithmic accuracy, the renormalization-scale dependence 
explicitly appears as a higher-order effect. 

Since we know the radiative factors only to NLL order, we use the Eqs. ( |32| , |36| , [39]) by replacing 
a s {k 2 ) (with k 2 = q 2 , (1 — z) 2 Q 2 , (1 — z)Q 2 ) with its NLO expansion in terms of a s (/x 2 ) and k 2 (cf. 
Appendix A), and we explicitly carry out the z and q 2 integrals by neglecting terms beyond NLL 
accuracy. We thus write the prompt-photons radiative factors as follows: 

Lou 2 ), %■ %) = exp {lniV ^(6b«.(A* 2 ) **N) + g^(b a s ^ 2 ) IniV, Q 2 /^ 2 ; Q 2 //4) 
V ^ V-fJ 



11 



+ 0(a s (a s \nN) k )} . 



(44) 



The functions g>W and resum the LL and NLL terms, respectively. These functions are 
different for the qq and qg partonic channels of Eqs. ( p6|) and (|27|), and are explicitly computed 
in Appendix A. We find: 



7$ (A) = (2C F - C A ) fc«(A) + C A ^(A/2) 



7 « 



(A) = C A /iW(A) + C F /iW(A/2), (45) 



and 



+ 



(2C F - C A ) /> (2) (A) + 2C a /i (2) (A/2) 



(46) 



+ 2Cf 7 Ca ln21n(l - 2A) + nb ° ln(l - A) - ^ Aln^ 

2ix% tcoq tx% Hp 



Cf 
nbn 



2A + ln(l -2A) 



+ 



2nbn 



21n(l - A) - ln(l - 2A) 



In 



Q 2 Q 2 



C A h^\\) + 2C F h { - 2 \\/2) 
° A In 2 ln(l - 2A) + 4C Ws ~ 3C F ln(1 _ A) 



(47) 



+ 



2nb 
Cf + Ca 
2irb n 



47r6 r 



CV + Ca , , Q 2 
1 Aln — 

7t6 up 



2A + ln(l - 2A) 



C h 



2nbn 



21n(l - A) - ln(l - 2A) 



1 2 

H 2 



where je = 0.5772 ... is the Euler number and bo, b\ are the first two coefficients of the QCD 
/3-function 



l\C A -AT R N f , nC A -10C A T R N f -6C F T R N f 



12vr ' h 24tt 2 

The auxiliary functions and in Eqs. and ( |4"E| , pT7|) are 

1 



fc«(A) 



27rfe n A 



2A + (1 -2A)ln(l -2A) 



/^(A) 



2A + ln(l-2A) + -ln 2 (l-2A) 



IE 

nb 



K 



111(1 " 2A) - 4^68 



(48) 



(49) 



2A + ln(l-2A) , (50) 



where K is the coefficient in Eq. (|35|). 

The results in Eqs. (fl"4]-|4~7D provide us with a theoretical description of soft-gluon resummation 
in prompt-photon hadroproduction at the same level of accuracy as for other hadropro duct ion 



processes, such as the production of Drell-Yan pairs ]3T], |32| or heavy quarks [^6], |2~8fl . These 
results can be used for detailed quantitative studies along the lines of Refs. [IT, In this paper 
we do not present numerical analyses and we limit ourselves to discuss the expected sign and size 
of the resummation effects. 

In the near-threshold region, radiation in the final state is kinematically inhibited. On physical 
basis, one thus expects that resummation of the ensuing logarithmically-enhanced corrections 
produces suppression of the cross section. This argument applies to hadronic cross sections, but it 
is not necessarely valid for partonic cross sections. The partonic cross section is what is left after 
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factorization of long-distance physics into the parton distributions. Since all-order resummation 
is in part automatically implemented in the definition of the parton densities, the remaining 
resummation effects can either enhance or deplete the partonic cross section. 



Among the various terms on the right-hand side of Eq. ([H]), some factors are smaller and 
some others are larger than unity. The exponent of the initial-state contribution in Eq.( 



is positive definite and, hence, >> 1 when N — > oo. The presence of this 'anti-Sudakov' 
form factor is typical of partonic cross sections that are evaluated by factorizing parton densities 
defined in the MS factorization scheme. In the case of the final-state contribution Jfj, no additional 
factorization has been performed. Therefore, when N — > oo the exponent in Eq. (|36|) is negative 
and J^r C 1 is a 'true' Sudakov form factor, as naively expected. The sign of the exponent in 
Eq. ( 



is not definite (D^ gry < 0, DjJ 



11^91 - "1" qg 

However, the contribution of A*" 1 ' 



dry • 



'(11 



Dqg'^q^ > 0) as expected for an interference term. 
N is subleading with respect to those of A a N and Jfj. 

From the overall inspection of the effect of the radiative-factor contributions to Eq. (|3T|), we 
infer that, in the case of prompt-photon production, the resummed partonic cross sections o^^L N 
and 7v i n Eqs. ( |26| ) and ( p7|) are both enhanced with respect to their LO approximations 

& q q-+gy Ni 0qg-*qy, n- Moreover, the enhancement in the qg partonic channel is larger than that in 
the qq channel. 

This conclusion can also be argued by a simplified treatment within the DL approximation. 
Inserting Eqs. (f4~2l , [43]) into Eq. fl3"ID , we obtain 



a 



(res) 

qg^-y, N 



= o. 



(0) 

qg^q-y, N 



a (res) 
a qq^y, N 



a (res) 
°W-»7, N 



~ a 



= a 



5$^, N exp | [2C F + 2C A - C F ] In 2 ivj 
exp|(C F + 2C A )^ln 2 iv} > 

expl[2C F + 2C F -C A }^ln 2 
exp|(4C F -C A )^ln 2 iv}> 



qg^qi, N > 



(o) 



a 



(res) 

qq—>~y, N 



11^11, N 



"(0) 
-(0) 

a qq^q~1, N 



.(0) 

qq— > <?7i N > 



exp\3(C A -C F )^-\n 2 N\ > 



a 



j.(0) 

qy^rh N 
(o) 

qq^q-y, N 



(51) 
(52) 

(53) 
(54) 

(55) 



The first, second and third terms in the square bracket on the right-hand side of Eqs. ([51], |53[) are 
respectively due to the initial-state factors A^, A^ and to the final-state factor Jfj. Note that, for 
a definite parton a, the initial-state enhancement Afj is larger than the final-state suppression Jjy 
(see the difference by a factor of two in the exponent of Eqs. ([42] , |43|) ). In the qg channel the final- 
state contribution In J% is thus overcompensated by In A q N and this leads to the enhancement in 
Eq. fl52"|) . In the qq channel, instead, it is the total initial-state contribution (In A^ + ln A 9 N ) that, 
owing to the colour-charge relation C F ~ C A /2, overcompensates In J 9 N . Finally, the enhancement 
in Eq. ([55]) is simply due the fact that the gluon colour charge C A is larger that the quark charge 



C F and, thus, A^ 



> A q N and J% 



> I 9 

> j n . 



Note that this conclusion directly applies only to the asymptotic limit iV — > oo or Ex — 
y/S/2. In the case of kinematic configurations of experimental interest, subleading effects and 
their dependence on the x-shape of the parton distributions and on the renormalization and 
factorization scale have to be carefully estimated. 
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4.3 The constant factors 



Expanding the resummed expressions in Eqs.([44]-|4~7|) to the first order in a s and using Eqs. (p6| 
ETft, we obtain 



a 



(res) 

qq— >7, N 



[a s (n );E T ,fi ,n F ,fi f ) = a a 



r (0) 



a s (/i 2 ) { 1 + 



a s (jj? 
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2C F -ho A )\nX 



+ [Ie^Cf - C A ) - (2C F - C A ) In 2 + nb - 2C F \n 



2E, 



T 



/4 



+ Ci;-l 7 (2E|//i 2 ;2^/^ 



InN 



(56) 



(7 



(res) 

<3ff^7> N 



M/ r)ji • nJ/,2) 



7T 



C F + C A ) In 2 JV 



qg^q~f, N 

+ ( 7i? (C , F + 2C A )-C J 4ln2 + ^C F -(C F + C A )ln 



3 „ , „ „ s , 2i?£ 

MniV 



+ (9(a 2 



/4 



(57) 



One can easily check that the logarithmic terms in these perturbative expansions agree with 
those that can be derived from the complete NLO analytic results of Refs. P], |42|, [|. From this 



comparison we can also extract the first-order constant coefficients CiiL,.- and Cj$ . We find 



99^7' 



{2C F -\ 



2 CUJ + 7 f [vr6 - (2C F - C A ) In 2j - - (2C F - C A ) In 2 



1 2 i cr 

+ -ir-^ + y(2C F --^) + -(2C F -C j4 )ln 2 2 
27 F C F - ^C F ) In % - vr6 In ^ , 

L fI F 



(5f 



where 



-C F -C A hi 2 



7|(^Cf + Ca) + 7f 
~ # g + |j (2C F + 19C A ) + \c F In 2 2 



(7f(Cf + C A ) - jC F - 7t6 ) In % - 7r6 In 



— (C7 F -2CU) In 2 



K n 



7 7T 2 



(59) 



(60) 



and the coefficient K is given in Eq. (|35|). 

The first-order coefficient and, indeed, all the perturbative coefficients of the constant 

(iV-independent) function C a b-+ 7 (a s ) in Eq. ( j30|) are produced by hard virtual contributions and by 
subdominant (non-logarithmic) soft corrections to the LO hard-scattering subprocesses. In both 
cases the structure of the external hard partons is the same as at LO. This justifies the all-order 
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factorization of C a b^(a s ) with respect to ^L^iV anc ^ ^o the radiative factor in the resummed 



partonic cross sections 

The inclusion of the iV-independent function C a b~ il {a s ) in the resummed formulae does not 
affect the shape of the cross section near threshold, but improves the soft-gluon resummation by 
fixing the overall (perturbative) normalization of the logarithmic radiative factor. 

We can explicitly show [28| the theoretical improvement that is obtained by combining the 
NLL radiative factor with the first-order coefficient Expanding the resummation formulae 

(26 , |27D in towers of logarithmic contributions as in Eq. ([I]), we have 



<rj cs) (« s ; E 2 , fi 2 , rf,) = a a s 1 + £ < [c„, 2 „ In 2 " N + (£|//4) In 2 "' 1 N 

I n=l 

+ Cn^Elj^Elltf) ln 2 "" 2 iV + 0(ln 2 "- 3 iV)]} , (61) 

where a s = a s (/i 2 ). The dominant and next-to-dominant coefficients c n . 2n and c n . 2n _i are con- 
trolled by evaluating the radiative factor to NLL accuracy. When the NLL radiative factor is 
supplemented with the coefficient we can correctly control also the coefficients c n>2n _ 2 . In 
particular, we can predict (see Appendix B) the large- N behaviour of the NNLO cross sections 
ogL, in Eq. (0) up to OQnN). 

Note also that coefficients c n 2n are scale independent and the coefficients c n 2n _i depend on the 



sole factorization scale //p. In the tower expansion (61), the first terms that explicitly depend on 
the renormalization scale /i (and on /ip, as well) are those controlled by c ni2n _ 2 . Their dependence 
on n is obtained by combining that of [E\j /ip, E\j /j 2 ) with that of the radiative factor at 
NLL order. The inclusion of the first-order constant coefficient thus (theoretically) stabilizes 
the resummed partonic cross section with respect to variations of the renormalization scale. 



Comparison with other processes: soft-gluon interfer- 
ences and QCD coherence 



Further insight on the underlying physics mechanism that leads to the resummed expressions (P6[ 
P7| ) can be obtained by comparing prompt-photon production with other hard-scattering processes. 

In the hadroproduction of a DY lepton pair (Fig. |l]a) of high mass Q 2 , the vicinity to the 
threshold region is measured by the inelasticity variable r = Q 2 /S, where v5 is the centre-of- 
mass energy. The Born-level partonic process that controls the cross section is qq annihilation. In 
vV-moment space, where the iV-moments are defined with respect to r, the Sudakov corrections to 
the gg-annihilation cross section are taken into account by a resummation formula analogous to 
Eq. (p6j ^7|) . Up to NLL accuracy, the corresponding radiative factor A dy ,n(Q 2 ) has the following 
explicit expression 

A DY , N (Q 2 ) = A%(Q 2 ) A%(Q 2 ) , (62) 

where A q N (Q 2 ) and A q N (Q 2 ) are the single-parton contributions^ given in Eq. (^). Each term 
A^r embodies multiple initial-state radiation of soft gluons, i.e. gluons that carry a small fraction 

*To simplify the notation, we drop the explicit dependence on a s and on the renormalization and factorization 
scale. Therefore, we use A a N (Q 2 ) = A^(a s (^ 2 ), Q 2 /^ 2 ; Q 2 //i|) and Jn(Q 2 ) = ^("sC" 2 ), Q 2 /^) throughout this 
Section. 
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Figure 1: Schematic representation of the structure of the soft-radiation factors in: (a) 
DY production, (b) DIS, (c) e + e~ annihilation, (d) prompt-photon photoproduciton and 
(e) prompt-photon hadroproduction. 



1 — 2 ~ l/N ~ (1— r) of the energy of the initial-state parton a. The factorized structure on 
the right-hand side of Eq. (|62l) implies that soft-gluon interferences between the two hard partons 
cancel to this logarithmic accuracy [44 . 

Note that this cancellation does not depend on the type of annihilating partons. In fact, when 
the DY pair is replaced by a colourless system, say, a Higgs boson, produced by gluon-gluon fusion, 



the resummed partonic cross section is controlled by a NLL radiative factor HQ, 45 



A Higgs , N (Q 2 ) = A 9 N (Q 2 ) A%(Q 2 



(63) 



which is again factorized in single-parton contributions. 

The presence of non-interfering Sudakov factors is typical of other processes dominated by 
hard scattering of two QCD partons, such as lepton-hadron DIS, e + e~ annihilation in two jets 
and prompt-photon photoproduction. 

In the case of inclusive DIS (Fig. |I|b), the hard-scattering scale Q 2 = —q 2 is given by the square 
of the space-like transferred momentum q and the relevant inelasticity variable is the Bjorken 
variable XBj = Q 2 /2Pi • q- The Born- level partonic process is lepton-quark scattering and, when 
the threshold region XBj — » 1 is approched, the corresponding radiative factor Adis, n{Q 2 ) in 



iV-moment space is [44, |46 



(64) 



The Sudakov factor A Q N (Q 2 ) is exactly the same as in the DY process. It embodies soft-gluon 
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radiation from the initial-state quark. Unlike in the DY process, however, in DIS the scattered 
initial-state quark fragments in the final state. Then the factor Jfj(Q 2 ) takes into account the 
fragmentation of the final-state quark into a jet of collinear and/or soft partons with a small 
invariant mass k 2 ~ Q 2 /N ~ (1 — XBj)Q 2 - The NLL expression of the jet mass distribution 
J%{Q 2 ) is given in Eq. 

Hadronic final states with two back-to-back jets produced in e + e~ annihilation at the centre-of- 
mass energy Q (Fig. |l|c) are also controlled by the jet mass distribution J^{Q 2 ) j43| . For instance, 



in the case of the distribution (l/cr) da/dT of the thrust T flT] , the Sudakov region is T — > 1. In 
this limit we have 1 — T ~ k\/Q 2 + k\jQ 2 , where k\ and k\ are the hadronic invariant masses 
in the two emispheres singled out by the plane orthogonal to the thrust axis. Considering the 
iV-moments A T ( e + e -^ N (Q 2 ) of the thrust distribution with respect to T, and taking the large- iV 
limit, one obtains [ETJ 

A T(e+e - )i7V (Q 2 ) = r N {Q 2 ) 4(Q 2 ) . (65) 

The factors J q and J 5 are the invariant-mass distributions of the two jets that originate from the 
fragmentation of the gg-pair produced by the e + e~-annihilation process at the Born level. 

The structure of the radiative factors in Eqs. ([62H651) easily explains the high-.&r behaviour of 
the prompt-photon cross section in photoproduction collisions (Fig. |l]d). This process, which can 
be regarded as a simplified case of the hadropro duct ion process considered throughout the paper, is 
discussed in Appendix C. In hadron-photon collisions the production of high- .E^ prompt photons is 
dominated at the Born level by the Compton-scattering subprocess q(xiP\) +7(^2) — * q(ps) +7(7?) ■ 
The all-order resummation of Sudakov effects leads to the radiative factor A q ^ qi in Eq. (|96]), 
whose NLL expression is given in Eq. ( |10U| ): 

A™(Q 2 ) = A%(Q 2 ) J q N (Q 2 ) ■ (66) 

The factor A q N (Q 2 ) takes into account sof-gluon radiation from the initial-state quark, while 
J%{Q 2 ) is the mass distribution of the jet produced by the collinear and/or soft fragmentation of 
the final-state quark. 

Note that high-E'y prompt-photon photoproduction can be regarded as a photon-hadron deep- 
inelastic scattering, where the space-like momentum transferred by the scattered photon is g M = 
P2 ~ ■ Since the high-i^ cross section is dominated by the kinematics configurations in which 
the prompt photon is produced in the central rapidity region, we have 2p2 ■ p — 1E\ and 2P\ ■ p ~ 
P\ ■ P2 = S/2. Thus, the hard scale is Q 2 = —q 2 = 2p 2 ■ p — 2E\ and the inelasticity variable 
analogous to the Bjorken variable is Q 2 /2Pi ■ q = Q 2 / (2P 1 ■ p 2 — 2P 1 ■ p) ~ AE^/S = x 2 T . Recalling 
that in Eq. ( |66"D we have Q 2 = 2E^ (see Eq. (|98])) and that iV is the moment index with respect to 
x\ (see Eq. (p5|)), we can thus straigtforwardly understand the complete analogy between Eq. 



and the expression (|64|) for the DIS radiative factor. 

The Sudakov corrections to prompt-photon hadroproduction (Fig. |I|e) are embodied in Eqs. ( pB| 
f27| ) through the radiative factor A c £^ dl (Q 2 ) . On the basis of the factorization of the right-hand 



side of Eqs. (|5"i2|-|6"6"D in terms of initial- and final-state single-parton contributions, one might expect 
that A c £~* d ' y can be obtained from the photoproduction result in Eq. (|66"D by simply including an 



additional initial-state factor A^. The NLL expression (0) for A^ d7 (Q 2 ) shows that this naive 
expectation is not correct. In fact, we have 



A a ^ d \Q 2 ) = A a N (Q 2 ) A b N (Q 2 ) 4(Q 2 ) A^ nt) " 7 (Q 2 ) . (67) 
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The presence of the NLL contribution A^ nt ' ) ab ^ dl Q n the right-hand side of Eq. fl67|) implies 
that the physical picture of the Sudakov radiative factors in terms of independent single-parton 
contributions is not valid, in general. As discussed in Sect. |2] and explicitly shown in Eqs. (|6^-p6|). 



this picture applies to processes dominated by hard-scattering of two sole partons, but it breaks 
down at NLL accuracy in the case of multiparton hard-scattering. The breakdown is due to 
interferences and colour correlations produced by soft gluons that are radiated at large angle with 
respect to the directions of the hard-parton momenta ||29|| . Soft-gluon interferences are present in 
the hard-scattering of three QCD partons as shown by Eq. (|67D, while colour correlations affect 
hard-scattering of more than three QCD partons . 



Owing to their large- angle origin, soft-gluon interferences are process dependent. In the case 
of prompt-photon hadroproduction they are taken into account by the factor /\(^ ab ^ d "f ^ w hose 
explicit NLL expression is given in Eqs. (]39| |4T|). 



Note that the coefficient D^_^ d in Eq. (|4l|) depends linearly on the colour charges of the hard 
partons, and the colour-charge dependence of A^ nt ^ ab ^ d ' y i s thus factorized at NLL accuracy. This 
suggests that the effect of the interference factor can be absorbed by a proper rescaling of the 
independent-emission factors A a , A b and J d . As a matter of fact, neglecting corrections beyond 
NLL order, one can check that the right-hand side of Eq. fl67|) can be rewritten as follows 



A^(Q 2 ) = A% /2 (Q 2 /2) A b N/2 (Q*/2) J d N/2 (Q 2 /2) . 

This equation has to be regarded as a manifestation of the colour-coherence properties of QCD 
emission |[j5|| . Soft gluons radiated at large angle destructively interfere. Their effect can thus be 
taken into account by Sudakov factors of independent emission in a restricted (angular) region of 
the phase space. 



6 Conclusion 

In this paper we presented the explicit expressions for the resummation of threshold-enhanced log- 
arithms in hadronic prompt-photon production, to next-to-leading accuracy. The simple colour 
structure of the diagrams contributing to prompt-photon production reflects itself in the simplic- 
ity of the resummed formulae. Fragmentation processes, furthermore, do not contribute to the 
Sudakov resummation at NLL level. In Mellin space, the resummed radiative factor factorizes in 
the product of three independent contributions for the initial and final coloured partons appearing 
in the Born process, times a simple factor describing the soft-gluon interferences between initial 
and final states. General coherence properties of large-angle soft-gluon radiation allow to further 
simplify the result: the interference contributions can be described, to the same degree of accuracy, 
by constraining the phase-space for independent emission from the coloured partons. The result- 
ing radiation factor can thus be written as the product of the three independent single-parton 
contributions, with a properly rescaled dependence on the Mellin-moment variable N. 

The formulae are given in terms of Mellin moments, and can be used for phenomenological 
applications by inverse-Mellin transforming to xt space. The problems related to this inversion 
are the same as those encountered in the resummation of the Drell-Yan or heavy-quark produc- 
tion cross-sections, and can therefore be solved with the same techniques [^TJ. All ingredients 
are therefore available for a phenomenological study of prompt-photon production including the 
evaluation of Sudakov effects with NLL accuracy. Such a study is in progress, and will be reported 
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soon. The calculations presented in this work, together with previous work on Drell-Yan and DIS, 
make it now possible to carry out global fits of parton densities with a uniform NLL accuracy in 
the large- x region. All of the processes that are used for these global fits, among which prompt- 
photon production plays a critical role, are now known theoretically at this level of accuracy. 



Acknowledgements We thank W. Vogelsang for useful discussions. 



Appendix A: NLL formulae for the radiative factors 



The logarithmic expansion of the radiative factors in Eqs. ([S^, |39|) can be computed as described 
in Refs. |32|, ^3|. The running coupling a s {k 2 ) with k 2 = q 2 , {l—z) 2 Q 2 , {l—z)Q 2 has to be expressed 
in terms of a s (/i 2 ) according to the NLO solution of the renormalization group equation: 



a s (k 2 ) 



a s (jj? 



1 + b a s ([i 2 ) \n(k 2 /fj, 2 



a s (n 2 



b 1 + b a s (fj, 2 ) \n(k 2 /fjL 2 



\n(l + b a s (^)\n(k 2 /^)) 



+ 0(a 2 s (fi 2 )(a s (^)\n(k 2 /ii 2 )r 



(69) 



where bo, b\ are the first two coefficients of the QCD /^-function, which are explicitly reported in 
Eq. (|4q) . Then the z integration can be performed with NLL accuracy by setting 



z"- 1 - 1 



-ea-z 



-7b 



/N) . 



Defining 
we find 

lnA a N (a s (ii 2 ),Q 2 / f i 2 ;Q 2 /i4) 
lnJ a N (a s (fi 2 ),Q 2 /^ 2 ) 

(int) ab— >dj / (,,2\ ^2/,,2\ 



A = 6 a s (/i 2 )lnA 



(70) 
(71) 



In A fc<?>(A) + h^(\,Q 2 /fi 2 -Q 2 /fi 2 F ) + O (a s {a s \nN) k ) ,(72) 
In AT /W(A) + /< 2 >(A, Q 2 /fi 2 ) + O (a s (a s lnN) k 



In A 



N 



M^),Q7// 



D (i) 

2-Kbn 



ln(l-2A) + C(a s (a s mA) fc ) 



(73) 
(74) 



where the LL and NLL functions h)p, ffp and hS 2 \ ff& are given in terms of the perturbative 
coefficients A^\A^\ B « in Eqs. $K 



h^(X,Q 2 /^,Q 2 /^ 2 F ) 



+ 



2nb X 
2nbl 
nb 



2A + (1-2A) ln(l-2A) 
2A + ln(l - 2A) + - ln 2 (l - 2A) 



(75) 



2vr 2 6 2 



ln(l - 2A) 
2A + ln(l - 2A) 
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+ 



Al 1 ) 



2tt6 



2A + ln(l-2A) 



(76) 

In— r-Aln — , 



/i 1} (A) 

/i 2) (A,g 2 //^ 2 ) 



Am 

a 

2vr6 A 

27lbl 



(1 - 2A) ln(l - 2A) - 2(1 - A) ln(l - A) 
ln(l - 2A) - 2 ln(l - A) + - ln 2 (l - 2A) - ln 2 (l - A) 



(77) 



+ ^Vln(l-A)-^V^ 



2n 2 b 2 



21n(l - A) - ln(l - 2A) 



ln(l - A) - ln(l - 2A) 



(7f 



+ 



2nbc 



21n(l - A) - ln(l - 2A) 



In 



Q 2 



2 ' 



Note that the functions /«(A) and /^(A, Q 2 //i 2 ) can also be written in terms of and 
as follows 



/1 1} (A) = h^\\/2)-h^{\) 



(79) 



/i 2) (A,Q 2 //^ 2 ) = 2 fff)(A/2, g 2 //. 2 ; 1) - h?(\, Q 2 /^ 2 ; 1) + ^ ± Ml ~ A) ■ (80) 



27r6 n 



Inserting the expressions (|72| , |73| , |74|) into Eq. (|31~1), and using the explicit form of the pertur- 
bative coefficients A^\A^\B^\ D { ^ dj in Eqs. fl3|, ||, ||D we obtain the results in Eqs. (@-|47|). 



Appendix B: NNLO partonic cross sections at large N 



According to the notation in Eq. (|61~1) , the large- N behaviour of the NLO cross section fr^* in 
Eq. ( |TT| ) is written as 



-(1) (t?2 2 2 2\ _ -(0) 

a ab-*-j,N\ nj T^ 5 fJ'FiH'f) — a ab-^d~i,N 



cg> In 2 N + 47 (E't/ii f ) In iV + ^ (E'^/fip, E%/tf) 



S ab )( Z?2 /,,2 



(aft) / j-,2 / ,,2 7-i2 /,,2\ 



+ 0(1/JV)_ , 

where the various coefficients can be read from Eqs. (|SB|, [57] ) 



Ci 9 ?^ — — ( 2C f Ca 

1,2 7T V 2 



51) 



52) 



1 

7T 
1 

7T 



7i?(4C F - C A ) - (2C F - C A ) In 2 + Tvb - 2C F In H^L 

Te (C f + 2C A ) - CU In 2 + - (C F + C A ) In ^ 

4 n F 



(83) 



c^{E 2 /^E 2 /^) 
c^iE'/^E*/^) 

C^ 1 are given in Eqs. (E3, |5J 



-C^(2E 2 T /^;2E 2 T /^ 
1 

7T 



C^{2El/^2E*/v 



F) ) 



54) 
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(2) 

Analogously, we can write the NNLO cross section <7„ b i^ 7 as follows: 



*(2) . , -2 2 2 2 i 



,(0) 

ab^d-y, N 
>6) 



c ( 2 f ln 4 iV + 47(^/^)ln a iV 



+ c y 2 a ° 2 \E T /^ F ,E^/^)\ii A N + 0{\nN) 



J5) 



The coefficients 02,4, 02,3, 02,2 can be calculated by expanding the resummation formulae fl26| , |27| ) 
to the second order in a s . We find 

2 



c 2,4 



1 

2 



c l,2 



(86) 



,(99) _(«?) 



4 9) (^//4) 



3tt 
2 



4 



(87) 



<8${E%l&,E%ll?) = \ [c^(E T /^ F )] 2 + c[f c^\e t /^E t /^) 



IT 



1e(±C f - |c A ) - (2C F - C A ) In 2 + + (2C F - 



+-b 

7T 



lE (^C F + 2C A ) -C A \n2 + ho F + (±C F + C A 



K 



1 

1,1 It. 



K , 2£f/ 
— m — — 



27rfe n 



where the coefficients 01,2, ci 5 i, c^o and if are given in Eqs. fl82| , |33| |J) and fl35|) . 

Our prediction for the coefficients in Eq. fl8"5] ) can be used to check future NNLO calculations 
of the prompt-photon production cross section. Alternatively, when these calculations become 
available, they can provide a highly non-trivial check of our NNL resummation. 
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Appendix C: photoproduction of prompt photons 



In hadron-photon collisions the inclusive production of a single prompt photon is due to the process 

H 1 (P 1 )+ 1 (P 2 )^ 1 (p)+X . (89) 

We use the same kinematics notation as in the hadropro duct ion case (cf. Sect. [|) and we write 
the prompt-photon photoproduction cross section integrated over 77 at fixed Ex as follows: 

da^ h \x T ,E T ) = t daW{x T ,E T ) \ + ( da^\x T) E T ) \ 

\ dEjr J jjacjjoiuc \ Et J p j n yijjg 

The hadronic contribution to the cross section is completely analogous to the right-hand side of 
Eq. (|10|) apart from replacing fb/H 2 ( x 2, Hf) with the parton distribution f b / 1 (x 2 , fi F ) of incoming 
photon. 

The second contribution on the right-hand side of Eq. (|90|) is due to point-like interactions of 
the incoming photon with high-momentum partons. The point-like cross section can in turn be 
decomposed in direct and fragmentation components 



■ dx I 5\x - j a a7 _ 7 (x,a s (/i 2 );E|,/i 2 ,/i|,/iy) (91) 



+ ^ j Q dZ ^ d ch( Z ^)) S ( X ~ ^^=^a7-»c(^,0! fl (/X 2 );E|-,/X 2 ,//|,,//J)| . 

The rescaled partonic cross sections a" a7 _> 7 and cr ai ^ c have perturbative QCD expansions similar 
to Eqs. (|TT|) and (12). In particular, for the point-like direct component we have 



2 ,.2 ,.2 ,.2^i _ a 2 



+ E «?(a* 2 ) ^ftM*! E t, A* 2 , /4, A*/) 

n=l 



(92) 



where the only non- vanishing terms at LO are those due to the Compton scattering subprocesses 

q + i->q + i, 9 + 7-*? + 7> (93) 

whose contribution to the cross section is 

*£L*r(*) = *gL#r(*) = t e{ ^=1= (l + ^ . (94) 

To perform soft-gluon resummation at high Et, we work as usual in iV-moment space by 
defining 

The resummation of the large- iV corrections to the iV-moments of the hadronic contribution 
in Eq. ( p0|) is exactly the same as for the hadroproduction case discussed in Sect. |j. Moreover, in 
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the large- N limit, the point-like contribution turns out to be dominant: the hadronic contribution 
involves the additional convolution with the photon parton density /f,/ 7 and this implies its sup- 
pression by a relative factor of 0(1/ N). We can thus limit ourselves to considering the point-like 
cross section. 

In the case of the point-like contribution, one can repeat the argument in Sect. ^4. 1| on the 
relative size of the fragmentation component and of the various direct subprocesses. Up to NLL 
accuracy, we then conclude that soft-gluon resummation in the photoproduction cross section (95) 



is controlled by the point-like direct channels q^y — > 7 and q^y — > 7. The all-order resummation 
formulae for the corresponding partonic cross sections are 

a^ l N (a s (v 2 y,E*,fi 2 ,fi 2 F ,ii 2 f ) = a 2 £jJL g7 , n C q7 ^ y (a s (n 2 ), Q 2 //j?; Q 2 /n%) 

■ a^v>*(/AQ 2 /V;Q7/4) , (96) 

<T^ 7jiV (a s (^ 2 );E|,/i 2 ,/i|,/iJ) = a^ 7j7V (a s (/i 2 );E|,,/i 2 ,/i^,/i5) , (97) 

where 



Q 2 = 2E 2 , (98) 

^=< TJ rmwr {7+5N) - (99) 



and n are the iV-moments with respect to x 2 of Eq. 



The radiative factor A q ^ qi and the iV-independent function C qi ^ in Eq. ( p7|) can directly be 



related to the analogous contributions A q ^ qi and Cg ff -+ 7 to the qg channel in the hadropro duct ion 



process. 



The radiative factor A^ - * 97 is obtained from the factorized expression ([H]) for A^^ 97 , namely 

99-*<n _ A? AS jQ /v(in1 
N — N N N N 



from A q ^ qi = A q N A 9 N J'^A^ t ^ q9 ~* q ' y , by switching off soft-gluon radiation from the incoming 



gluon. This amounts to set Ca = in A 9 N and ^jy 1 ff->97 . Using the explicit formulae in 
Eqs. ( P2"| and ( |3"9"| , [Hp , this implies that up to NLL accuracy we can neglect both Afy and 



Ajy 9S_>?7 and we have the simple result: 

A^> s (/i 2 ), Q 2 //, 2 ; Q 2 //4) = A q N (a s (fi 2 ), Q 2 / ^ Q 2 //40 J q N (a s (^), Q 2 /V) • (100) 

Note that no soft-gluon interference factor A^ nt ' ) appears in Eq. ( |100|) . Prompt-photon pho- 
toproduction at threshold is dominated by an underlying hard-scattering that involves only two 
hard partons and then, in agreement with the general discussion in Sect. |2], soft-gluon interferences 
have to cancel. 

The explicit NLL expansion of Eq. ( |1UU[ ) gives 
Q 2 Q 2 ' 



AIT 97 « s (^ 2 ), \; \) = exp {lniV g^(b Q a s (^ 2 ) IniV) + g g 2 )(b a s (fi 2 ) IniV, Q 2 //i 2 ; Q 2 / ^ 



+ 0(a s (a s lnN) k )} ; (101) 



where the LL and NLL terms and g^ 2 ' are expressed in terms of the auxiliary functions 
and of Eqs. © and (0): 

9 <»{\) = C F hV{\/2), (102) 
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9$ i A 



g 2 g r 



97 l '" /i 2 >F, 



2 C F fc< 2 > (A/2) + 4 ^^ -3C f ln(1 _ A) 



C 



g 2 c* 



7TOo /Up 7T0, 



47r6 n 



A + ln(l-A) 



In 



/i 2 



(103) 



The ^-independent function C qi ^(a s ) has the following perturbative expansion 



(/Ag 2 /V;g7/4) = ^ 7 l 7 (g 2 //4)+E (^) n ^t 7 (g 2 /^ 2 ;g 2 /^) • 

(104) 

Note that the first-order coefficient Cjff does not depend on the renormalization scale. Its 
explicit expression is obtained from that of C*ggL 7 by setting Ca = and bo = in Eq. (|59|) : 

^WV^)^,{i 7 I + 5 7E - ±ln2-i| + |! + iV2-( 7£ -?)l I1 |} . (105) 



Expanding the resummation formula fl96|) in powers of a s we can derive the large- N behaviour 
of the NLO and NNLO cross sections and of Eq. 

At NLO we find 



-(1) (Tp2 2 2 2\ _ ~(0) 



where 



In 2 iV + c?J>(Ey^) In iV + eg? + 0(1/JV) 



lii)(Tr^ /,,2 



(97)/ 7712 / ,,2 



(106) 



c! 9 2 7) = ^-C F , ^(^/^) - ic F ( 7 s + 7 - In ^ ) , c§\Ey&) = l -C^{2El/f F ) 
2tt it \ 4 n F J TV ' 

(107) 

This result agrees with the large- N limit of the NLO analytic expressions computed in Ref. [[5J. 
At NNLO we predict 



,(o) 



S<rr) (j?2 /,,2 \ i„3 



+ c? 2 7) £|//i 2 ) In 2 AT + C(ln A/") 



where 



c 2,4 



„(97) 
c l,2 



8tt 2 



_ / 3 , 2£f 



71 r 



(108) 

(109) 
(110) 



J,Ti)( E T E ¥ 
2 ' 2 1^' 7", 



2 
2tt 2 



C^l^E^/fil) + (C7 F + tt&o) (7s + 7 - In 



2E% 



K 

1 

-K 



In 



2£ 2 



(111) 



and the coefficient K is given in Eqs. (|3q) . 
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